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The recent paper claims that mean characteristics of cheiotic orbits differ from the corresponding 
values averaged over the set of unstable periodic orbits, embedded in the chaotic attractor. We 
demonstrate that the alleged discrepancy is an artifact of the improper averaging: Since the natural 
measure is non-uniformly distributed over the attractor, different periodic orbits make different 
contributions into the time averages. As soon as the corresponding weights are accounted for, the 
discrepancy disappears. 

PACS numbers: 05.45-a 



Recent Rapid Communication ^ compares properties 
of periodic orbits (UPOs) embedded into a chaotic at- 
tractor to those of chaotic trajectories on this attractor. 
Analysis is based on numerical data, and culminates in 
the statement: "time-averaged properties along a set of 
UPOs and a set of chaotic orbits with finite lengths are 
totally different from each other" . It is further conjec- 
tured that "the time averages of the dynamical quanti- 
ties along UPOs with the same period of the Poincare 
map have a limiting distribution with nonzero variance" . 
In this comment we show that under the proper aver- 
aging procedure, there seems to be neither the "total" 
difference between the averages, nor an argument for the 
above conjecture. 

Due to ergodicity, the value of time average for an ob- 
servable A converges to ^ = J A{x)fi{x)dx where /i is 
the natural measure, and integration is performed over 
the whole attractor. Approximation of A by summation 
over the set of UPOs - apparently the method employed 
in [H - provides correct results only in exceptional cases 
when fi is uniformly distrubuted over the attractor: for 
linear mappings like the Bernoulli map, symmetric tent 
map etc. In general, however, the density of the natural 
measure varies along the attractor. As a consequence, a 
chaotic orbit does not walk uniformly over the attract- 
ing set: it visits certain regions relatively often or stays 
there relatively long. Contribution of such regions into 
the time averages is larger than of those visited seldom. 
As shown in [Ij (see also Chapter 9.5 of the textbook 0]), 
non-uniform distribution of the natural measure can be 
recovered from the properties of UPOs embedded in the 
attractor. In particular, for invertible two-dimensional 
maps (and hence for three-dimensional flows which in- 
duce such maps), the weight with which an UPO con- 
tributes to the time average, is inversely proportional to 
the largest eigenvalue of the corresponding fixed point. 
Below we demonstrate that taken the weights into ac- 
count, the discrepancy between the mean values from the 
chaotic timeseries and the mean values from the UPOs 
disappears. 



We restrict ourselves to the first example considered 
in [l]: attractor in the celebrated Lorenz equations 
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with CT = 10, r = 28,6 = 8/3 [1, and take the same 
observable: the value of the variable z. In fact, already 
in [5!| fractal characteristics of the Lorenz attractor were 
evaluated with the help of the properly weighted UPO 
data. Our results are obtained by averaging z for each 
orbit from the complete set of 111014 UPOs with TV < 20 
turns in the phase space around either of the attractor 
"wings" . For each iV < 20 we also computed 10^ seg- 
ments of chaotic trajectories with N turns, and calcu- 
lated values of {z) for every such segment. 

Comparison of eigenvalues for different UPOs with the 
same N discloses strong inhomogeneity in the distribu- 
tion of natural measure. Already among 186 orbits with 
TV = 11 the largest (4618.57) and the smallest (415.59) 
eigenvalues differ by the factor of 11.1. At = 19 there 
are 27594 UPOs, and the ratio between the extremal 
eigenvalues is 77.2. Accordingly, the contribution of the 
"most unstable" UPO is hardly discernible, compared to 
the contribution of the "least unstable" one. 

As seen in the left panel of Figlll taking the weights 
into account shifts and re-shapes the distribution of mean 
values. The solid curve shows the bell-shaped distribu- 
tion of (z) for segments of chaotic orbits with length 20. 
The dashed curve shows the histogram obtained by sum- 
mation of (z) from all UPOs of the same length. Simi- 
larly to Fig.l of [21], the maxima of these two curves are 
shifted with respect to each other. The expectation val- 
ues of (z) for these two distributions are distinctly differ- 
ent: (z) equals 23.555 for chaotic trajectories and 23.420 
for summation over UPOs. This difference, however, al- 
most vanishes for the histogram which incorporates the 
weights of UPOs: the dotted curve in FiglTJa) is much 
closer to the solid curve, and has (z)=23.554. In the 
case of shorter orbits with length A^ = 11 reported in [ij 
the same effect takes place: ensemble of chaotic orbits 
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FIG. 1: Time averages of z{t) for chaotic trajectories and 
UPOs of the Lorenz equations. 

(a) Probability density for orbits of length Ai'=20. Solid 
line: average values from 10^ segments of chaotic trajectories. 
Dotted line: weighted means over all UPOs with A'^ = 20. 
Dashed line: summation (without weights) over all UPOs 
with N = 20. 

(b) Standard deviation of {2:) for periodic (pluses) and chaotic 
(crosses) orbits with A'^ turns. Dashed lines: power-law fits. 

predicts (z) =23.562, which is definitely distinct from the 
value 23.420 obtained by summation over UPOs, but is 
much closer to the value 23.550 yielded by proper sum- 
mation with weights. 

Obviously, the distribution obtained from finite seg- 
ments of chaotic orbits cannot be identical to the distri- 
bution produced by the set of UPOs of the same length. 
A'' rotations in a piece of a chaotic orbit do not necessar- 
ily constitute a visit into the neighborhood of an UPO 
of the same length: they can consist of several visits into 
vicinities of the shorter UPOs, be a part of the passage 
near the longer UPO, etc. Therefore - at least in the 
checked range of iV - distributions based on chaotic or- 
bits are broader than their counterparts built from the 
UPOs data. As shown in the right panel of FiglU the 
breadth of both distributions decays as a power of 
For chaotic segments it approaches the law 3.04 iV~°-^^ 
(which nearly coincides with the values from whereas 
the standard deviation for a distribution from the set of 
UPOs with properly assigned weights decays as 
Of course, the employed values of N are, at best, moder- 
ate, and one cannot judge on the asymptotical properties 
of the dependence. At any rate, however, these data un- 



ambiguously show that in the range iV < 20 there are no 
arguments against eventual convergence to the 5-shaped 
distribution [7|. 

Since for small and moderate values of N the distribu- 
tions are typically broad, it hardly makes sense to discuss 
whether an attractor of a particular set of equations is 
accurately modeled by a single UPO of the given length: 
some of the UPOs are definitely inappropriate, whereas 
some others may prove to be good. Average values com- 
puted along the UPO match the averages along chaotic 
orbits either in a pathological case when all UPOs except 
one possess giant eigenvalues (and, hence, negligible sta- 
tistical weights), or if one deliberately chooses the UPO 
whose characteristics are close to that of the whole en- 
semble. For the latter, however, the ensemble (or, at 
least, its representative parts) should be examined, so 
there is hardly a gain in the computational efficiency. 

The last remark concerns estimates of the topological 
entropy in Jj. It is known (see e.g. Q) that at r = 28 
a chaotic orbit can make not more than 25 consecutive 
turns around one wing of the Lorenz attractor before a 
jump to another wing. Once each turn in the half-space 
X > is coded by " 1" and each code in the complemen- 
tary half-space a; < is coded by "0", all binary strings 
with length < 25 are met in the code of a sufficiently 
long chaotic trajectory. The only missing periodic or- 
bits are those whose symbolic labels consist exclusively 
of ones or of zeroes. Accordingly, the number of UPOs 
with the length N for < 25 is given by the recur- 
sive formula K{N) = (2^ - 2 - J2j j K{j)j /N, sum- 
mation being taken along all divisors j of N. How- 
ever, estimate of the topological entropy as htop = 
limjv^oo sup A^~^ log K{N) = log 2 is applicable only for 
the newborn attractor at r = 24.06 ... in which all sym- 
bolic strings of arbitrarily large N are encountered. This 
is not the case for r—28: upwards from N=25 the tree of 
symbolic sequences is "pruned" , and the growth of num- 
ber of UPOs as a function of N may become slower. In 
any case, the range of orbit length < 14 employed 
in (ij for the evaluation of the htop is far too short and 
hardly appropriate for reliable estimates. 
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